Exploring exchange mechanisms with a cold atom gas 
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Fermionic atoms trapped in a double well potential are an ideal setting to study fundamental 
exchange mechanisms. We use exact diagonalization and complementary analytic calculations to 
demonstrate that two trapped fermions deliver a minimal model of the direct exchange mechanism. 
This is an ideal quantum simulator of the Heisenberg antiferromagnet, exposes the competition 
between covalent and ionic bonding, and can create, manipulate, and detect quantum entangle- 
ment. Three trapped atoms form a faithful simulator of the double exchange mechanism that is the 
fundamental building block behind many Heisenberg ferromagnets. 

PACS numbers: 67.85.Lm, 03.65. Ge, 03.65.Xp 



Recent experimental advances allow investigators to 
confine up to ten atoms in a trap and address their quan- 
tum state [TJ 12]. This innovation enables the Heidel- 
berg group to isolate two distinguishable fermions in a 
one-dimensional well, and tune the interaction strength 
to induce fermionization [5], presenting a unique oppor- 
tunity to study the fundamental physics of short range 
repulsion [4lfTT|. This could allow experimentalists to 
realize an analog to the Stoner model for itinerant fer- 
romagnetism [T3[T3]. However, many real- life solids are 
best described by spins that are localized in real-space. 
These are commonly described by the Heisenberg model 
that predicts the direct and double exchange mechanisms 
behind the emergence of antiferromagnetism and ferro- 
magnetism respectively. Here we take advantage of the 
ability of the Heidelberg group to trap either two or three 
fermions in the double well potential to realize the first 
exact quantum simulators for the direct exchange and 
double exchange mechanisms. This allows us to not only 
study the minimal building block of Heisenberg mag- 
netism, but moreover build a quantum simulator to ex- 
pose the competition between covalent and ionic bond- 
ing, and study quantum entanglement manipulation and 
detection. 

The direct exchange mechanism has previously been 
realized in double quantum dot systems [ T4HT7] , and cold 
atom gases in an array of double well potentials [TBHIO] ■ 
We now exploit the experimental flexibility of the Hei- 
delberg group to isolate just two atoms and expose the 
phase behavior. We use exact diagonalization to deliver 
the full energy spectrum, and a perturbative approach 
to gain an intuitive understanding of the phase behavior. 
The flexibility of the setup allows us to build the first 
quantum model of the fundamental covalent and ionic 
bonding mechanisms in molecules and crystals, allowing 
us to address the long-standing question of their relative 
contributions to chemical bonds [211 . A thorough under- 
standing of the exchange energy also allows us to define 
a new protocol to create, control, and detect quantum 
entanglement P^I23| . 

The isolated double exchange mechanism has not yet 
been realized experimentally. A faithful quantum simu- 
lator of the double exchange mechanism is not only im- 



portant for understanding Heisenberg ferromagnets but 
also describes the 90° superexchange mechanism, and can 
be extended to larger system through a quantum virial 
or cluster expansion [24] [25j. Through exact diagonal- 
ization and a complementary perturbative analysis we 
demonstrate how the cold atom gas will probe the dou- 
ble exchange mechanism with changing barrier height, 
interaction strength, and the ellipticity of the external 
trapping potential. This modifies the degeneracy of the 
ground state leading to a quantum phase transition that 
we expose with a statistical tunneling probe. 



I. FORMALISM 

We start from two-component Fermionic atoms with 
each species indexed by a pseudospin a G {tii}- The 
atoms are trapped by an external harmonic potential 
with Hamiltonian 7? = -V"^ /2 + [w\{x'^ + y'^)+uj'^^z'^]/2 + 

Vb exp(-a;B2:^)+g(r^-r;)n^(r^)n;(r;), setting h = m^ 
1 throughout. The parabolic trapping potential has a 
variable ellipticity that we associate with a length scale 



= iA 



The Gaussian barrier VBexp{—i^BZ ) will 



split the system into two wells. Throughout this paper we 
set its width with wb = Swy . We use two complementary 
techniques to analyze the system: firstly exact diagonal- 
ization to expose the full energy spectrum, and secondly 
perturbation theory to deliver an intuitive description. 

In exact diagonalization |27| , we work in the eigenbasis 
of the non-interacting Hamiltonian. We express the or- 
bitals in this basis as linear combinations of the Gaussian- 
type orbitals (GTOs) [28] of the harmonic trapping po- 
tential without the central barrier. Labeling these GTOs 
by the standard quantum numbers {rix, "-y, '^z}, we in- 
clude < (rix, %) < 4 and < n^ < 50. Including basis 
functions with higher n^ is necessary to capture the effect 
of the central barrier on the orbitals. 

To calculate the orbitals of the non-interacting system 
including the central barrier we first evaluate the non- 
interacting Hamiltonian matrix in the GTO basis. The 
matrix elements of the central Gaussian barrier can be 
conveniently expressed in terms of hyper-geometric func- 



tions f5S1. Having calculated the Hamiltonian matrix, 
we diagonalize it to obtain a list of orbitals to use in 
the subsequent calculation that incorporates the effect of 
interactions. 

We construct the 10,000 Slater determinants with low- 
est non-interacting energy to use as the many-body basis 
for our exact diagonalization calculation. To include in- 
teractions we evaluate the four-center integrals in the ba- 
sis of our orbitals and construct the Hamiltonian matrix 
using the Slater-Condon rules f^T" . Finally, we diagonal- 
ize the matrix to find the ground and excited states of 
the system. 

In a cold atom gas the Feshbach resonance is used 
to tune a contact interaction strength V{r) from the 
repulsive through to the attractive regime. To gener- 
ate a positive scattering length for a potential with a 
small effective range we use the square well potential 
g{r) = —g9{R — |r|) with a radius R = 0.4a||. This is 
significantly less than the width of the central barrier 
2/y/LJB ~ 0.9aj| so that atoms localized in adjacent wells 
will not interact, and moreover we ensured that the en- 
ergy of the states converged in the limit i? — > 0. The 
well depth g was tuned to generate a positive or nega- 
tive scattering length a = R[l — ta.n{Ry/g) / R^], and 
was constrained to < Ry/g < 4.49 to confine at most 
one bound state. However, the inclusion of the bound 
state leads to many molecular states crossing the low- 
est energy open channel (see Fig. flla)). This requires 
us to adiabatically track states between calculations per- 
formed at neighboring scattering lengths a^, a^+i, that 
can be uniquely followed by connecting states with the 
correct spin, inversion symmetry, angular momentum L^, 
and the greatest wave function overlap ('!/'m(fli)IV'"(^J+i)) 
where m, n are state indices. 



II. DIRECT EXCHANGE 

We first study two fermionic atoms trapped within 
the double well potential. This is a minimal quantum 
simulator of the direct exchange mechanism, exposes the 
competition between covalent and ionic bonding, and al- 
lows us to probe quantum entanglement. These three 
applications demand a thorough understanding of the ex- 
change energy that we study with two complementary ap- 
proaches: exact diagonalization and perturbation theory. 
The exchange energy can be probed directly in experi- 
ment by applying a magnetic field gradient and measur- 
ing the tunneling rate of atoms from the system ^ |3] . 
The general interacting system with no central barrier 
has been studied by Busch [6] so we first address a per- 
turbative barrier before focusing on a high barrier. 

Shallow barrier: We analyze the system analytically 
with first order perturbation theory. For the 5 = 
state we start with two opposite spin atoms in the Gaus- 
sian non-interacting ground state of the trapping poten- 
tial with no central barrier (f>^(r), and for 5 = 1 with 
two equal spin atoms with one excited into a state with 
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(c) Qubit control protocol 
(i) Prepare qubit in state |4t) 



(ii) Before evolution: qubit in state |4,t) 

(iii) Evolve with interactions: singlet-triplet oscillation 



(iv) After t = TT/E: SWAP 
Qubit in state |t4) 
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(v) Tilt to resolve and measure 




(vi) After t = 2iT/E: identity 
Qubit in state |4,t) 

±/\i 

(vii) Tilt to resolve and measure 




FIG. 1. (Color online) (a) Exchange energy from exact diag- 
onalization (black lines), with the open channel in red, lowest 
S = 1 state in blue, and perturbation theory (green dashed 
line), (b) Lower: Normalized exchange energy of the S = G 
level as a function of relative well depth e from exact diago- 
nalization (blue crosses) and from the Hubbard Hamiltonian 
(red line). Upper: The density profile (green) in the trap 
(magenta) at three values of the tilt and interaction strength, 
(c) Schematic of the experimental protocol to study quantum 
entanglement. 



a node either in the longitudinal or transverse direc- 
tions depending on whether un or uj±_ is smaller. Next 
we introduce the central barrier that within first order 
perturbation theory induces a change in the energy of 
^^ Jdr|(/)(7(r)pVBexp(— cjbz^), and scattering that ac- 
cording to mean-field theory increases the energy by 
47ra/ dr|(/).|-(r)p|0^(r)p. This yields an exchange energy 
(energy difference between the 5 = and 5=1 open 
channel states) of —A -I- oa/a\\ with A = min(a;||,a;_L) — 
Vb{['^||/('^|| +wb)]^/^ - [w||/(w|| +u^yf'/'^}, and overlap 

matrix element o = ^J 2 / ttu ±\\ ~ A^/2Vb / T^iijJw +wb)]- fu- 
creasing both the barrier height and interaction strength 
lowers the exchange energy. However, exact diagonaliza- 
tion confirms that the exchange energy is negative for 
any positive scattering length, so that a negative scatter- 
ing length is required for a 5 = 1 ferromagnetic ground 
state in conformance with the Lieb-Mattis theorem [29]. 
High harrier: We now focus on a high barrier that 



presents a unique opportunity to study a minimal Hamil- 
tonian of localized spins. We initially exploit exact diag- 
onalization to deliver the exact ground and excited states 
before developing an analytic theory to gain a heuristic 
understanding of the system. Finally, we demonstrate 
the versatility of the setup to study the competition be- 
tween covalent and ionic bonding, and quantum entan- 
glement. 

Results for a barrier with Vb = 4w|| arc shown in 
Fig. [Ha). The exchange energy flips sign across unitarity, 
so that the red S* = level is the ground state for a > 0, 
and the blue S — 1 level for a < 0. On approaching 
unitarity from positive scattering lengths many molecu- 
lar bands anti-cross the open channel from < a|| /a < 3, 
removing the energy gap to the other states. Raising 
the central barrier increases this energy gap, thereby re- 
ducing the number of molecular band crossings. This 
should facilitate experiments that adiabatically transit 
across unitarity. Although the two atom system cannot 
fall into a bound state without violating energy conser- 
vation, the closure of the energy gap precludes using this 
region to explore quantum entanglement. However, the 
Super- Tonks regime is free from the molecular bands, 
making it the ideal venue to use the system as a quan- 
tum simulator and study quantum entanglement. 

We next develop a complementary analytical expres- 
sions for the ground state energy to gain an intuitive 
understanding of the system's behavior. We model 
the system by the Hubbard Hamiltonian t J^a-i^ha'^R'^ + 
cLcl<t) + eE^(cLcL^ - 4<TCRa) + g(c[^c[^CL^CLt + 
cj^^c]^, cr^cr^), where c^^ creates a particle localized in 

the left-hand well, cj^^ in the right-hand well and e de- 
notes the relative depth of the wells. Diagonalization in 
the high-barrier limit t <C (e, go) yields an exchange en- 
ergy — (1 -I- e^an/16a^o^)t^a||/4ao, where we use WKB 
perturbation theory [30] to derive the parameters t = 
32cj|| exp(-2v/7rVB/wB)/7r2, and o = 3uj^^/{2tt)^^^. 
We apply a magnetic field gradient dB/dz adding a term 
AH = az to the Hamiltonian, where a — gj/iediJ/dz 
and gj is the g-factor. This gives a relative well depth 
of e ~ [log(2a;BVB/w^)/a;B]^/^a. We find that the per- 
turbative analysis matches well to the results of exact 
diagonalization around unitarity, as shown in Fig. [l[a,b) . 

Covalent versus ionic bonding: The covalent bonding 
in a crystal is driven by the negative exchange energy of 
electrons being shared between neighboring atoms. Con- 
versely, in ionic bonding the constituent ions carry dif- 
ferent electronegativities that drives a displacement of 
electron density resulting in bond polarization and con- 
comitant bonding. The cold atoms present an ideal sim- 
ulator of the chemical bonding mechanisms, where the 
external potential portrays the atomic pseudopotential, 
the fermionic atoms represent the valence electrons, and 
the tilt e the relative electronegativity tuning from e = 
(covalent) to |e| ~ wy (ionic character). Fig. [iTb) shows 
that the exchange energy is —(1 -|-e^an/16a^o^i^a||/4ao. 



where in the absence of the tilt the first term — t^a||/4ao 
corresponds to the covalent bonding, and the second term 
further lowers the energy in the presence of the tilt so cor- 
responds to the ionic contribution. On introducing the 
tilt Fig. [ljb:i-ii) demonstrates how the bond becomes 
polarized, with the increasing density in the lower well 
^ a'^,t^e/32a'^o^ lowering the net energy quadratically 
^ e^ conforming with Pauling's definition of electroneg- 



ativity [m [3T] . The lowering of the energy and increas- 
ing polarization of the bond could be detected by tun- 
neling atoms out of the two sides of the trap. Finally, 
Fig. [ljb:iii) shows that with a < the atomic density 
would be counter-intuitively pulled into the higher well. 
Quantum entanglement: As the experimental setup of- 
fers a unique level of experimental control we are well 
positioned to explore quantum entanglement. The qubit 
control protocol is shown in Fig.[lVc). We define the first 
qubit state as the up-spin localized in the left-hand well 
and the down-spin in the right-hand well, denoted ||t)j 
Fig. [IJc:ii & c:vi), and the second qubit state conversely, 
Iti), Fig. [T|c:vi). In Fig. [T|c:i) the system is prepared 
into a qubit basis state in the non-interacting regime in 
the presence of a magnetic field gradient gj/ied-B/dz > 
t/au in the Zeeman regime so that the two spins are 
driven in opposite directions [2 . With the configuration 
in Fig. [TJc:i) we form the qubit state Hf) in Fig. [l|^c:n). 
When the tilt is removed the qubit state ||t) is a su- 
perposition of the singlet and triplet states. With in- 
teractions the system will evolve under Rabi oscillations 
between the singlet and triplet states in Fig. [IVc:iii) with 
a Rabi period 2Tr/E, with E = —t^a\\/Aao. A duration of 
tt/E corresponds to a SWAP rotation 32J into the qubit 
state Iti) shown in Fig. [ljc:iv), and a duration 2tt/E 
an identity rotation into the qubit state |4,t) shown in 
Fig. 1 ^c:vi). To detect the qubit states in Fig. [ljc:v) and 
Fig. 1 |c:vii) we apply a strong magnetic field gradient in 
the Paschen-Back regime to empty the right-hand well 
but leave the left well occupied since it is shielded by 
the central barrier. The escaped atom carries a definite 
spin, resolving the system onto the qubit basis, and the 
remaining atomic spin can then be measured separately 
following the protocol in Ref. [5] to identify the qubit 
state. Finally, we note that di?/de|g=o = and therefore 
the Rabi period is insensitive to the controlling tilt pa- 
rameter |33j , making this an ideal opportunity to explore 
entanglement in a flexible, clean, and stable environment. 



III. DOUBLE EXCHANGE 

Two trapped atoms with a positive scattering length 
always yields a 5 = ground state ,29], so to realize 
the simplest possible ferromagnetic ground state we turn 
to a three-atom system. For a sufficiently strong central 
barrier this delivers a quantum simulator of the double 
exchange mechanism for ferromagnctism: two atoms are 
localized in the lowest level of each well, and a third itin- 
erant atom in a higher energy level couples their spins. 
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TABLE L Lowest energy states for the three atom system and their energy calculated in perturbation theory. Both the 5=1/2 
and S — 3/2 states are shown, and also the available longitudinal and transverse modes. 
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FIG. 2. (Color online) The phase diagrams predicted by ex- 
act diagonalization (left) and perturbation theory (right) at 
three different trap ratios. The phase corresponding to the 
each color is identified in each diagram, with "L" indicating 
a longitudinal phase and "T" a transverse phase. The ar- 
rows show the critical interaction strength predicted by the 
Hubbard model. 



With three atoms there are two possibilities for the high- 
est occupied orbital: either with a node longitudinally 
(riz = 1), or two degenerate modes with a transverse 
node (rix.y = !)■ To study the emergence of magnetic 
correlations we use both exact diagonalization and the 
perturbation theory. We adopt the same perturbation 
theory as introduced in the direct exchange section for 
a small central barrier and weak scattering, except now 
for the three-atom states shown in Table lU To orient the 
discussion we focus initially on the realization of itiner- 
ant ferromagnetism before introducing a central barrier 
to consider the double exchange mechanism. 

Exact diagonalization delivers the phase diagrams in 
Fig. [2] that compare well with the perturbation the- 
ory predictions. At Vb = the system enters the 
S = 3/2 state at a/a\\ « 0.45, which compares fa- 
vorably to the perturbation calculation estimate that 



a/an — -\/2/7r/3 ~ 0.84 and a prediction that we ex- 
trapolated from the data of Liu of a/a\\ ss 1.1 [8 . We 
note that the exact theory has a lower critical interaction 
strength than the perturbation theory prediction, which 
reflects the situation of the itinerant case [351 [3S] ■ The 
critical interaction strength is predicted to be indepen- 
dent of trap ellipticity uj±/tjj\\ by both exact diagonal- 
ization and the perturbation theory calculations. In the 
high barrier potential limit we follow the example from 
the direct exchange calculation and model the system 
with the Hubbard Hamiltonian, predicting a crossover at 
a/a^\ — •\/7r/2/3 « 0.42. This is in direct agreement with 
the results from exact diagonalization. 

When W|| < u!± it is energetically favorable to oc- 
cupy longitudinal modes and when wn > uj± the trans- 
verse modes are preferred. However, the introduction 
of a central barrier favors longitudinal states (n^ — 1) 
with a node across the barrier, with the crossover at 
Vb = ■\/T+"wb7w|[(1+W||/wb)(wjj —uj±) for both S = 1/2 
(at a — 0) and S = 3/2. A weak central barrier 
forces the atoms apart reducing the effective interac- 
tion strength, meaning that the boundary between the 
5 — 1/2 and S — 3/2 phases has a positive slope of 
y^{uj\\ + uj±) /24y^. With a larger central barrier the 
longitudinal S = 1/2 state and transverse S = 3/2 state 
share a phase boundary with a negative slope giving rise 
to a characteristic notch. 

The phases are distinguished not only by their spin 
quantum number but also the degeneracy associated with 
the orientation of the node in the highest occupied or- 
bital. The changing degeneracy means that the phases 
with the same spin quantum number cannot evolve into 
one another so there is a quantum phase transition be- 
tween the two. In a cold atom gas both total spin and 
symmetry are conserved, so to probe the phase diagram 
a tunneling measurement is proposed f2j, with the gas 
starting from four atoms, tilting the trap, and tunneling 
down to the three atom configuration with the lowest en- 
ergy. At the critical interaction strength the S ~ 1/2 and 
S = 3/2 states will be formed with equal likelihood, but 
taking account of the degeneracy will boost the forma- 
tion of the doubly degenerate states, presenting an ideal 
tool to characterize the phases. 

The high barrier perturbation theory provided a per- 
fect description for the phase behavior predicted by exact 



diagonalization theory in that Umit. This demonstrates 
that the ferromagnetism is driven by double exchange, 
consistent with the observed exponential decay of the in- 
teraction matrix element between one atom localized in 
the left and another in the right-hand well. However, the 
high barrier perturbation theory is unable to describe the 
system with a low central barrier, where the requirement 
to have orbitals with a longitudinal node is incompatible 
with the node-less "transverse" wave ftmctions found in 
exact diagonalization and shown in Table |lj which are of 
itinerant nature. Therefore, the crossover from "trans- 
verse" to "longitudinal" wavefunctions in the phase dia- 
gram denotes the underlying ferromagnetic correlations 
changing from itinerant to localized double exchange. 



between covalent and ionic bonding. The experimental 
flexibility of the cold atom gas also presents the ideal 
arena to study quantum entanglement. Trapping three 
fermions delivers the first faithful realization of the dou- 
ble exchange mechanism. This can provide insights into 
Heisenberg ferromagnets and 90° superexchange, and be 
built up to larger lattices through a cluster expansion. 

The exchange mechanisms presented here give a tanta- 
lizing insight into the broad range of effects that can be 
explored in the double well potential. The inclusion of a 
third potential well would allow investigators to study 
the Kramers-Anderson superexchange mechanism |34| 
behind many antiferromagnets. Trapping more atoms 
should reveal an even/odd effect of flipping between an- 
tiferromagnetic and ferromagnetic ground states. 



IV. DISCUSSION 

The double well system is the perfect playground to 
study the exchange mechanisms of strongly interacting 
fermions. We have studied the direct exchange mecha- 
nism that is the fundamental building block of Heisen- 
berg antiferromagnetism, and exposes the competition 
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